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The impedance and the radiation field in the surrounding hot plasma are calculated for 
an electrically short antenna that consists of two spherical conductors excited through thin 
wires in phase opposition. In the calculation the pressure tensor is replaced by a scalar 
pressure. A discontinuous model of the ion sheath is used. 

The losses due to the radiation of electromagnetic and electrostatic waves are calculated 
and are expressed in terms of equivalent series resistances. The operation of resonance 
probes is discussed. It is shown that their resonant frequency is well below the electron 
plasma frequency if the probe radius is much larger than the Debye length. The significance 
of this result to both past and future ionospheric rocket probe experiments is pointed out. 
The limitations of the present treatment are discussed. 



1. Introduction 

Antenna impedance measurements in the iono- 
sphere during rocket flights show the presence of a 
much larger resistive component [Whale, 1963] than 
is predicted by the type of antenna theory in which 
the surrounding plasma is treated as cold. In 
principle the interaction of an antenna with a hot 
and very tenuous plasma should be treated by solving 
a combination of the collisionless Boltzmann equa- 
tion with Maxwell's equations. The mathematical 
difficulties of solving these equations and satisfying 
the appropriate boundary conditions are, however, 
formidable, and therefore the problem is considerably 
simplified in this paper by the use of certain approxi- 
mations. The first of these is the use of a scalar 
pressure [Spitzer, 1962, p. 24] and the resulting 
neglect of Landau damping [Landau, 1946]. The 
problem is further simplified by the rather crude 
additional assumption that the plasma surrounding 
the antenna is sharply bounded and that outside 
the boundary the plasma is uniform. The effects 
of the earth's magnetic field and dissipative effects 
are neglected. For other applications of the same 
type of approximations the reader is referred to 
papers by Gould [1959], Fejer [1964], and Nickel, 
Parker, and Gould [1963]. 

These approximations are applied here to a simple 
dipole antenna that consists of two spheres excited 
in phase opposition through wires whose capacity is 
neglected. Two types of waves can be excited by 
such a dipole antenna in the uniform plasma: elec- 
tromagnetic waves and acoustic-type, longitudinal 
electrostatic waves. The distance between the 
centers of the two spheres is assumed to be much 
shorter than the wavelength of electromagnetic 
waves but much longer than the wavelength of elec- 
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trostatic waves and than the diameter of the spheres. 
The term "electrostatic waves" is used in this paper 
to describe the plasma oscillations discussed by 
Tonks and Langmuir [1929] and further analyzed 
by Landau [1946]. 

The present analysis and its conclusions differ 
from those of Whale [19G3], who also considered 
the excitation of electrostatic waves by introducing 
the concept of an isotropic pressure. Whale assumed 
the existence of interaction between the fluctuating 
quasi-electrostatic field (the near field of the electro- 
magnetic waves) and the electrostatic waves through- 
out the uniform plasma. In the present treatment 
the only interaction is taken to occur at the inner 
boundary (in reality in the sheath). Within the 
uniform plasma, in the absence of a magnetic field, 
the electrostatic wave and the electromagnetic wave 
(or in the present limit the quasi-electrostatic field) 
propagate independently and therefore cannot inter- 
act in terms of the hydrodynamic approximation. 

Since the completion of the present work a number 
of closely related papers [Cohen, 1961 ; Cohen, 1962a; 
Cohen, 1962b; Hessel and Shmoys, 1962; Chen, 
1963; Harp, 1963; Pavkovich, 1964; Wait, 1964] 
came to the author's attention. Of these the work 
of Cohen, Hessel, and Shmoys, and Chen are re- 
stricted to the calculation of the electromagnetic 
and electrostatic wave fields of point or line current 
sources in hot uniform plasmas. Wait considers a 
finite slotted-sphere antenna immersed in a hot 
plasma and his work is more closely related to the 
present one. Pavkovich, in contrast to all other 
workers, does not restrict himself to the concept of 
an isotropic pressure and gives numerical results of a 
calculation based on the collisionless Boltzmann 
equation, for the radiofrequency properties of a one 
dimensional plasma sheath. His results were modi- 
fied approximately for spherical geometry by Harp. 
This modification is closely related to the present 
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work and reaches similar conclusions about the 
theory of the resonant probe used by Miyazaki et 
al. [I960]. The relationship between the present 
work and the work of others will be further clarified 
in the body of the paper. 

2. Excitation of Electrostatic Waves 

The antenna described in the introduction radiates 
both electrostatic and electromagnetic waves. In 
the present section only one of the two spheres of 
the antenna will be considered at first and the near 
field of the electromagnetic wave field will be approxi- 
mated by the quasi-electrostatic field calculated 
from Poisson's equation. The impedance of the 
whole antenna, excluding the electromagnetic radia- 
tion resistance, is then simply taken to be equal to 
twice the impedance of a single sphere. 

The quasi-electrostatic field calculated in this sec- 
tion for a distance large compared to the distance 
between the spheres but small compared to the 
electromagnetic wavelength will be extrapolated in 
the following section to much greater distances. In 
this manner the field of the electromagnetic wave and 
the electromagnetic radiation resistance of the an- 
tenna, neglected in this section, are calculated in the 
following section. 

Let R be the radius of one of the spheres and r the 
distance of an arbitrary point from the center of the 
sphere. It is assumed that the unperturbed electron 
concentration vanishes for r<^R } and has the con- 
stant value N for r^>R. This is admittedly a rather 
artificial model of the unperturbed ion sheath. It 
effectively assumes an abrupt potential barrier (this 
could be visualized as a hypothetical double layer 
formed from infinitely heavy positive and negative 
ions) at r=R which prevents the penetration of 
electrons inside the sphere at r=R. It is clear that 
the radial component of the mean electron velocity 
must vanish at r=i£. Immediately inside the sheath 
there is assumed to be a spherical conductor whose 
(quasi-electrostatic) perturbation potential is taken 
to be a harmonic function of the time. (In principle 
the radius of the conductor could be taken as smaller 
than the radius of the spherical plasma-free cavity 
without essentially modifying the analysis; this will 
not be done here.) Since there is no mean radial 
motion of the electrons at r=R, there is no fluctuat- 
ing charged surface layer there and therefore both 
the perturbation potential V and its normal deriva- 
tive must be continuous at r—R. Inside the plasma 
the equation 



mNbv/bt=NevV-yKTvn 



(1) 



is taken to be valid where m is the mass, e is the 
absolute value of the charge, N is the unperturbed 
number density and n is the perturbation in the 
number density of electrons (factors e iat are taken 
for granted in n and in V), K is Boltzmann's con- 
stant, T is the temperature, and where the ratio y of 
the specific heats is taken as 3 [Spitzer, 1962]. 
A combination of the divergence of (1) with the 



equation of continuity (satisfied by the velocity v of 
the electron gas) 



leads to 
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A combination of (3) and Poisson's equation in 
spherical coordinates and with spherical symmetry, 
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where a 2 =(u N * — w 2 )/u 2 , «^=(eW/e ra) 1/2 is the elec- 
tron plasma frequency and u=(yKT/m) 1/2 is the 
velocity of electrostatic waves in the high frequency 
limit (co>>coat). 

The general solution, which vanishes at infinity, 
of the differential equations (4) and (5) has the form 
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where that part of the solution which is associated 
with the constant Ci, describes the electrostatic 
wave, whereas the part associated with C 2 describes 
a quasi-electrostatic field which is a good approxima- 
tion to the electromagnetic wave field within dis- 
tances very much shorter than a wavelength from 
the antenna. 

The boundary condition v=0 at r=R yields with 
the aid of (1) 
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Substitution of n and V from (6) and (7) into (8) 
leads to the relation 



a— - 
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Substitution of <7 2 from (9) into (6) and (7) leads to 
the following relationship between the perturbation 
in the number density and the perturbation potential 
V at r=R: 



V(R) 
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(10) 



Similarly, the effective, complex capacity of the 
conductor, defined here as the ratio of the charge 
on the conductor to the potential V(R) is given by 



a, 



V(R) 
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Equations (6), (7), (9), (10), and (11) represent 
the main results of this paper, and a discussion of 
their significance follows here. 

Equation (11) may also be used to calculate the 
impedance Z ett =2(io)C eU )~ 1 of the antenna consisting 
of both spheres. The energy loss due to the radia- 
tion of electromagnetic waves has not as yet been 
included in the analysis; this will be done in a sub- 
sequent part of the present paper. The resistive 
part of the impedance Z cff represents only the energy 
loss due to the "radiation" of electrostatic waves. 

It is convenient to express the impedance Z eff as 
a function of the parameters F=2(47reoi?co) _1 (twice 
the free space reactance of a sphere of radius i?), 
8=Ru~ l u N (the ratio of the radius R to 7 1/2 times the 
Debye length), and \p=u/u N (the ratio of the fre- 
quency to the electron plasma frequency). The 
result is 



1-8^(1 _^)-'/ 



0' 



It is clear from (12) that Z eff is purely imaginary 
for \p<^\ or ox^cotv. The absence of a loss term is 
explained in this case by the absence of propagating 
electrostatic waves and the neglect of Landau damp- 
ing. It may be seen from (12) that the impedance 
becomes infinitely high at the plasma frequency and 
vanishes at the frequency where 1 — 8\// 2 (l — \l/ 2 )~ i/2 =0. 
More will be said about the significance of this 
frequency later. At very low frequencies, Z ett 
= — %Y(l + <5) -1 , and thus the impedance is smaller than 
the free space impedance by a factor (1 -f-5)" 1 . It is 
understandable that the effective capacity of the 
sphere is larger at low frequencies than the free 
space capacity because the alternating electric field 
does not extend to infinity, but is confined to within 
the plasma sheath, which in the case of large 8 is 
much smaller than the radius. Figure 1 shows, for 
i^<0, the ratio of the imaginary part of Z eff (the real 
part is zero for iA<l) to the magnitude Y of the free 
space reactance as a function of the ratio \//=w/oo N - of 
the frequency to the electron plasma frequency for 
different values of the ratio 5 of the radius to 7 1/2 
Debye lengths. The curve for 8= oo represents the 
well-known cold plasma approximation in which 
thermal motions are neglected and the plasma is 
treated as a dielectric with an effective dielectric 
constant e (l— co^ 2 /^ 2 )- 

For ^>1 it is convenient to write (12) in the form 



rsy-i)- i/2 . rq+sy 2 ) 
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(13) 



Figure 1 also shows the real and the imaginary parts 
of Z ef f for ^>1 (in units of Y) as functions of i^, for 
different values of 8. The case of 8= oo again rep- 
resents the approximation in which thermal motions 
are neglected. The curves of figure 1 show the 
presence of a resistive component of the impedance 
for finite values of 8. The resistive, ohmical com- 
ponent represents the energy loss caused by the 
radiation of electrostatic waves. 



Im (Zetf/Y) 

Re(Z eff /Y) 




Figure 1. The ratios of the real (interrupted line) and imagi- 
nary (solid line) parts of the impedance to the magnitude of 
the free space reactance as functions of the ratio, &, of the 
radio frequency to the electron plasma frequency. 

3. Radiation of Electromagnetic Waves by 
Two Spherical Conductors Excited in 
Antiphase 

The analysis of this paper can be extended to in- 
clude radiation losses due to electromagnetic waves. 
The analysis of electromagnetic radiation by spheres 
oscillating in antiphase is, of course, somewhat artifi- 
cial because in practice antennas resemble cylindrical 
conductors. The analysis of spheres oscillating in 
antiphase has, however, the advantage of simplicity 
and it illustrates the nature of the problem rather 
well. It would be relatively simple to extend the 
present analysis to cylindrical conductors, at least in 
an approximate manner. An exact formal solution 
for a slotted sphere antenna has been given by Wait 
[1964]. 

Since the distance D between the two spheres is 
taken to be much smaller than the electromagnetic 
wavelength (and at the same time /2<<Z7), the 
quasi-electrostatic field at a distance r^>^>D on the 
line connecting the spheres is 2C 2 D/r*, where C 2 /r is 
that part of the potential given by (7) which could 
be regarded as an approximation to the radiation 
field of a spherical radiator at distances much smaller 
than the electromagnetic wavelength. At distances 
much longer than a wavelength and in the plane 
that perpendicularly bisects the line connecting the 
two spheres, the magnitude E of the radiation electric 
field of the dipole antenna is then given by 

E= C 2 Dk 2 /r= (OJD/r) [ (co 2 - a N *) j c 2 ] . (14) 

The magnitude H of the radiation magnetic field 
is given by 

H=(e C 2 D/r)[(u 2 -co N *y' 2 /ca>] (15) 

where k is the wave number and c is the velocity of 
electromagnetic waves in vacuum. Equations (14) 
and (15) were derived by fitting the quasi-electro- 
static dipole field to the radiation field of a dipole 
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[Stratton, 1941 J in a medium with a dielectric con- 
stant € (1 — co^/co 2 ). Equations (6), (7), and (9) can 
then be used to express C 2 in terms of V(R), the 
alternating potential of one of the spheres, or in 
terms of the current supplied to the spherical con- 
ductors I=i(aC et iV (R) where (7 eff is given by (11). 
The Poynting vector EXH may then be expressed 
in terms of /, and the total radiated power P may 
be found by integration over a very large sphere. 
The result is 



P={<OTre,)- l D 2 uC-\a> 2 -u N *) l/2 P 



(16) 



and therefore the resistive component of the antenna 
impedance (the radiation resistance) due to the 
radiation of electromagnetic waves is 



P^-(67T6o)- 1 C 2 COC- 3 ( W 2 -CO^)l/2. 



(17) 



The resistive component due to electrostatic waves 
is given by (13) as 



.=r 



<3QA 2 -1)- 1/2 
1+W-l)' 



(18) 



Equations (17) and (18) show that p em is zero at the 
electron plasma frequency and increases rapidly 
with frequency whereas p es is infinitely high at the 
plasma frequency and decreases rapidly with fre- 
quency. 

It is interesting to note that the expression (17) 
for p em is the same as the expression for the radiation 
resistance of a short antenna of length D (whose 
capacity is entirely at its extremities) situated in a 
medium with a dielectric constant € (1 — u N 2 /a) 2 ). 
This is a significant result in view of the fact that 
in the calculation of the series reactive part of the 
antenna impedance the plasma cannot be replaced 
by a medium with a dielectric constant e (l — w N 2 /a) 2 ). 

It is interesting to compare the ratio p em \Pes 
resulting from (17) and (18) with (73) of Wait for 
the corresponding ratio of the power P v radiated in 
the form of electrostatic waves to the power P e 
radiated in the form of electromagnetic waves by 
a slotted sphere antenna. Since Wait's equation 
(73) presupposes k v a^>^>\ in Wait's notation, the 
corresponding assumption k v R^>^>l permits the 
neglect of unity in comparison with h 2 (\p 2 — 1) in the 
denominator on the right of the present equation 
(18). Equations (17) and (18) of this paper then 
result in the equation: 



Pes _ 

Pem 



0>^2 u 1 



c D 2 R 2 kj 



(19) 



where Wait's definition of k e is used. The factor 
2a~ 4 in Wait's equation (73) has been replaced by a 
factor 3D~ 2 R~ 2 in (19) but otherwise the two equa- 
tions are identical. The similarity of the present 
result to that of Wait has thus been demonstrated 
for the special case k p R^>y>l. 



4. Application to Resonance Probes 

The results of the previous section may be used to 
draw some tentative but very important conclusions 
about the behavior of resonance probes [Miyazaki 
et al., I960]. In such probes the change, caused by 
the application of a radiofrequency voltage in the 
direct current collected by a Langmuir probe, is 
measured as a function of the radiofrequency. It 
has usually been accepted [Miyazaki et al., 1960] as 
an experimental fact that the change in the collected 
current shows a resonant increase at the plasma 
frequency. 

In this section the point of view is taken that the 
change in the collected direct current is due to 
rectification caused by the nonlinear characteristic 
of the Langmuir probe. The amplitude of the radio- 
frequency variations in the collected current will be 
proportional to the fluctuations n(R) in the number 
density (in a more accurate treatment the fluctuation 
in temperature would also have to be taken into 
account but these would in any case be proportional 
to the fluctuations in number density) and therefore 
(10) should roughly describe the characteristics of 
a resonance probe. Using the parameters 5 and \p> 
(10) may be written in the form 



n(B)=V(B)^{ 



l 



5xf> 2 



~{\-W\ 



i)" 1 - (20) 



Equation (20) shows that n(R) becomes infinitely 
large at the frequency for which 



^={2- 1 r 2 [-l + (l+45 2 ) 



2\ 1/21 11/2 



(21) 



This shows that resonance occurs near the electron 
plasma frequency OA^l) only when 5«1> for 
large values of d the resonance occurs well below the 
plasma frequency. Figure 2 shows (eu 2 /o) N 2e Q )\n(R)/ 
VCR) I as a function of \p for two values of 5 and 
illustrates that resonance occurs well below the 




Figure 2. The dimension-less quantity (eu 2 /co N 2 e ) |n(R)/ 
V(R)|, which is proportional to the number density fluctua- 
tions, n(R), produced by an rf voltage of given amplitude V(R), 
as a function of the ratio, \F, of the radio frequency to the 
electron plasma frequency . 
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plasma frequency when the probe radius is much 
larger than the Debye length. The resonant fre- 
quency given by (21) is the same as the previously 
discussed frequency where the impedance of the 
probe vanishes. 

The present analysis clearly leads to the conclusion 
that the resonant response of the probe does not 
occur at the plasma frequency, as laboratory and 
space experiments are claimed to show, but always 
below the plasma frequency. Large errors could 
thus occur in rocket investigations of electron con- 
centration in which the resonant frequency was 
identified with the plasma frequency, especially with 
probes whose radii were much larger than the Debye 
length. Fortunately a relatively small sphere (R = 
1 cm) was used in the ionospheric experiments [Aono 
et al., 1962] reported so far. 

A discussion of the physical nature of this resonance 
sheds some light on the reasons for the occurrence 
of the resonance well below the plasma frequency. 
It is clear from (7) that the electric field of the probe 
consists of two parts. One is simply a quasi- 
electrostatic field, which at short distances is the 
approximate form of the radiation field for fre- 
quencies above the plasma frequency whereas the 
other is the field associated with an electrostatic 
wave which is evanescent at frequencies below the 
plasma frequency. 

At very low frequencies the quasi-electrostatic 
field becomes very small compared to the field of t he 
evanescent electrostatic wave. This means that 
the alternating charge on the conducting sphere is 
almost perfectly shielded by a suitable (continuous) 
modification of the sheath. As the frequency is 
increased, the shielding becomes less perfect and 
a potential C 2 jr appears outside the sheath. The 
present theory shows that the outside field opposes in 
phase the field within the sheath. With increasing 

! frequency a situation is reached eventually where the 
potential drop outside the sheath just balances the 
potential drop inside so that no exciting voltage is 
required on the conductor and resonant oscillations 
occur. 

This explanation of the resonance is strongly 

I supported by the full agreement, in one special case, 
of the present approximate calculation with Landau's 
[1946] more accurate analysis utilizing the collision- 
less Boltzmann equation. Landau considered the 
penetration of an alternating electric field, that is 
normal to a plane boundary, from empty space on 
the one side of the boundary into a uniform plasma 
on the other side. He showed that at frequencies 
slightly below the plasma frequency the field asymp- 
totically approaches a value which is equal to 
(I-coa^/ to 2 ) -1 times its value at the boundary. 
The field deep inside the plasma is thus out of phase 
with the field at the boundary. Equations (6), (7), 
and (9) of the present paper Lead to the same asymp- 
totic value if R and r-R both approach infinity and 
at the same time satisfy the inequality i2^>^> r-R. 
Landau also showed that the asymptotic value is 
approached according to a simple exponential law 
expressed by his equation (45). Essentially the 



same exponential law with the same exponent (- a r 
in equation 6) results from (6), (7), and (9) of this 
paper. 

The above conclusions about the behavior of the 
resonance probe described by Miyazaki et al. 
[1960] are in qualitative agreement with most of the 
conclusions of Harp [1963] who used an approximate 
modification of results, obtained with the aid of the 
collisionless Boltzmann equation by Pavkovich 
[1964], to spherical geometry. The most important 
difference between the present treatment and that of 
Harp [1963] and Pavkovich [1964] is the neglect of 
Landau damping by the present treatment; this 
damping appears to reduce the resonance effect 
drastically for probes whose radius is smaller than 
about twice the Debye length [Harp, 1963]. 

5. Conclusions 

The simple analysis of this paper leads to certain 
interesting and significant results about the behavior 
of antennas in plasmas and about the interpretation 
of observations with the aid of resonance probes in 
plasmas. It is shown that resonance does not take 
place at the plasma frequency and that previous 
measurements made by the resonant probe method 
may have to be reinterpreted. In principle both 
the electron concentration and the electron tem- 
perature could be determined by simultaneous 
measurements of the resonant frequencies of two 
resonance probes of different size. Alternatively the 
impedance of a single probe could be measured at 
two different frequencies. A theory that is more 
accurate than that of the present paper or than 
other presently available theories, would, however, 
be required for the accurate interpretation of such 
measurements. 
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